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Abstract 

Peakons are singular, soliton-like solutions to nonlinear wave equa- 
tions whose dynamics can be studied using ordinary differential equations 
(ODEs). The Degasperis-Procesi equation (DP) is an important exam- 
ple of an integrable PDE exhibiting wave breaking in the peakon sector 
thus affording an interpretation of wave breaking as a mechanical colli- 
sion of particles. In this paper we set up a general formalism in which 
to study collisions of DP peakons and apply it, as an illustration, to a 
detailed study of three colliding peakons. It is shown that peakons can 
collide only in pairs, no triple collisions are allowed and at the collision 
a shockpeakon is created. We also show that the initial configuration of 
peakon-antipeakon pairs is nontrivially correlated with the spectral prop- 
erties of an accompanying non-selfadjoint boundary value problem. In 
particular if peakons or antipeakons are bunched up on one side relative 
to a remaining antipeakon or peakon then the spectrum is real and simple. 
Even though the spectrum is in general complex the existence of a global 
solution in either time direction dynamics is shown to imply the reality of 
the spectrum of the boundary value problem. 

1 Introduction 



The prototypical example of PDEs admitting peaked solitons is the family 

Ut - U^xt + (h+ l)uUx = bUxU^x + UU^xx, (1-1) 

often written as 

rrit + rrixU + bmUx = 0, m = u — u^x, (1-2) 
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which was introduced by Degasperis, Holm and Hone [T], and shown to be 
Hamiltonian for ah values of 6 [2]. The most studied cases are the Camassa- 
Holm (CH) equation {b = 2), and the Degasperis-Procesi (DP) equation [Tl[3] 
{b = 3). For 6 > these are the only values of b for which the equation is 
integrable, according to a variety of integrability tests [3HS]- The case & = is 
relevant for another reason; this case provides a regularization of the inviscid 
Burgers equation that is Hamiltonian and has classical solutions globally in 
time [?]■ The 6- family admits an appealing geometric interpretation as Euler- 
Arnold equations on the space of densities m{x)dx^ for the group of orientation- 
preserving difFeormorphisms DifF(S'^) [5]. 

In order to discuss peakon solutions one needs to develop the concept of 
weak solutions. Because of the role that the Lax pair formalism plays in the 
theory we will define weak solutions in such a way that the PDE in question is 
the compatibility condition of weak Lax pairs. This prompts the x member of 
the Lax pair to be viewed as an ODE with distribution coefficients while the t 
equation of the Lax pair is viewed as a (isospectral) deformation of the former. 
We subsequently need to rewrite the PDE itself as a distribution equation. 
To this end we observe that the formulation (|1.2p suffers from the problem 
that the product mUx is ill-defined already in the case of continuous, piecewise 
smooth (in x) u{x,t), since the quantity m — u — Uxx is a measure with a 
non-empty singular support at the points of non-smoothness. To make matters 
worse that measure is in addition multiplied by the function Ux which has jump 
discontinuities exactly at those points. This problem can be resolved easily; one 
instead rewrites as 

(1 - dl)ut + {b+l-dl) dx (i u") + dx ul) = 0. (1.3) 

The case 6 = 3 is of particular interest to us. Then the term is absent 
from equation (|1.3p and in that particular case one requires only that u{-,t) G 
L[^Q^(R); this means that the DP equation can admit solutions u that are not 
continuous [9HTT]. 

Multipeakons are weak solutions of the form 

n 

u(a:,0 (1-4) 

i=l 

formed through superposition of n peakons (peaked solitons of the shape e"'^'). 
This ansatz satisfies the PDE (jl.3|) if and only if the positions (xi, . . . , a;„) and 
momenta (mi, . . . , rn„) of the peakons obey the following system of 2n ODEs: 

n n 

ifc = X!"^*^"'"^""'^''' rizfc = (6-l)TOfc^mi sgn(a;fc-Xi)e~l^''~^'l. (1.5) 
1=1 1=1 

Here, sgna; denotes the signum function, which is +1, —1 or depending on 
whether x is positive, negative or zero. In shorthand notation, with ^/(x)) 
denoting the average of the left and right limits, 

(/(x)) = i(/(x-) + /(x+)), (1.6) 
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the ODEs can be written as 

Xk = u{xk), mfc = -(5- l)TOfc(w2,(xfe)). (1.7) 

In the CH case (6 = 2) this is a canonical Hamiltonian system generated hy h = 
^ rrij m/c e~l^j~^'=l, for which xjs and nijS are canonical positions and 

momenta. In the DP case (b = 3) this is a non-canonical Hamiltonian system 
with the Hamiltonian H — '^j ^^'^ ^ non-canonical Poisson structure 

given in [12]. 

It is important to distinguish the case of pure peakons (initial r7ij(0) > 0) or 
pure anti-peakons (initial mj{0) < 0) from a general case of multipeakons (no 
restriction on the signs of mj(0)). Pure peakons have peaks, pure anti-peakons 
have troughs while multipeakons contain both peaks and troughs. 

The relevance of multipeakon solutions is that they provide a concrete model 
for wave breaking [13l[T4]. For more information on the wave breaking phe- 
nomenon for this class of wave equations the reader is referred to [?] and [TMlQ] . 
In the CH case the distinction between pure peakons or anti-peakons and multi- 
peakons does not result in a serious departure from the inverse spectral formulas 
for pure peakons. Indeed, explicit formulas for the n-peakon solution of the CH 
equation were derived by Beals, Sattinger and Szmigielski |2Qj and then ex- 
tended to n- multipeakons in jl4j|2T] using inverse spectral methods and the 
theory of orthogonal polynomials. The situation for the DP equation is con- 
siderably different. The analysis of pure peakon solutions for the DP equation 
was accomplished by Lundmark and Szmigielski [35J[53] using inverse spectral 
methods and M.G. Krein's theory of oscillatory kernels |24]. In short, in these 
papers, it was shown that in the DP case, when working with pure peakon or 
pure anti-peakon solutions, the concept of total positivity plays a fundamental 
role, for example, implying that the spectrum involved is positive and simple. 
For this reason going beyond the pure peakon sector of the DP will not be 
as straightforward as in the CH case which remains self-adjoint in the whole 
multipeakon sector. The DP spectral problem, by contrast, is manifestly non- 
selfadjoint. Yet, in addition to a general interest in modelling the wave breaking 
mechanism, there is another reason for studying multipeakon solutions of the DP 
equation: in |11] Lundmark introduced a new type of solution, a shockpeakon 
solution, which he showed for the case n = 2 gives a unique entropy weak solu- 
tion originating from the peakon-antipeakon solution. Therefore, multipeakon 
solutions can also provide us with an additional insight into the onset of shocks. 

The paper is organized as follows. In Section[2l we set up the formalism for an 
arbitrary number n of multipeakons and elaborate on general forms of peakons 
with index 1 and n. We emphasize the role of the boundary value problem 
and its adjoint, both associated with the a;-member of the Lax operator. In 
Section[3l we undertake a detailed study of three multipeakons. We establish the 
analytic character of solutions in Lemma 13.111 and describe the main properties 
of colliding pairs, culminating in Theorem 13.171 describing the creation of a 
shockpeakon at the collision. Section S] is devoted to analysis of the spectrum 
of the boundary value problem, in particular we prove a signature-type Lemma 
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14.11 relating the signs of masses mj (0) of colliding peakons to the real parts of 
eigenvalues. Finally, in Section [5] we classify different asymptotic in t behaviour 
of three multipeakons in terms of the sign configurations of the initial masses 
m,(0). 



2 Lax pair and the multipeakon spectral prob- 
lem 

It was shown in [Tj that the DP equation admits the Lax pair: 

(a, - a,,,)^' = zm*, ^t^[z'\l-dl)+U:,-ud,]-9. (2.1) 
In particular if u is given by the multipeakon ansatz (jl.4p . 



m 



(2.2) 



and equations (|1.7p for 6 = 3 follow from the (distributional) compatibility 
of equations (|2.ip . The boundary conditions consistent with the asymptotic 
behaviour of '5 read: 



as X ^ — oo, 



is bounded as 



(2.3) 



To see how the implementation of these conditions leads to an isospectral prob- 
lem we will trace back the most important steps in analysis of Lax pair for 
peakons. For more details the reader is referred to ^23j. We start in the region 
X < Xi lying outside of the support of the discrete measure to. There, the 
first equation in the Lax pair can easily be solved and the boundary condition 
implemented by ^'(x) = e^. When Xk < x < Xk+i, we have 

*(x) = Afc(z)e" + Bfc(z)+Cfc(z)e-^ 1 < fc < ti. (2.4) 

The coefhcients Ak{z), Bk{z), Ck{z) are polynomials of degree k in z given by 




fc]-) Vie/ 



p-1 



e 3 3+1 




where ('^p"') is the set of all p-element subsets I = {ii < 
For a; > ccn we will drop the subscript n, thus 



(2.5) 

< ip} of {1, ...,n}. 



■^{x) = A(z)e^ + B{z) + C(z)e^^ 



(2.6) 
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In the case of interest for us, namely n — 3, the coefficients A, B, C can be 
written exphcitely 

A{z) = 1 - [mi + TO2 + ms] z 

+ [mim2(l - e"=i-"=2)2 + TO2m3(f - e"=^"^^)2 + miTO3(f - e"^!-"^^)^] 

- [mim2m3(l - e^i-="^)2(l - e'=^-'''f]z^ (2.7a) 
i3(z) = 2 [mie^^^ + m2e'^^ + mae^^] z 

- 2[TOim2(l - e^i"="^)2e^^ + m2m3(f ~ e'=^-''^fe'=^+mim3{l - e^^-''--'fe'-']z 
+ 2[TOiTO2m3(l - 6^1-^^)2(1 - e^^-^^)2e^3]z3, (2.7b) 
C(z) = ~ [mie^^i + 77126^^^ + mae^^^'] 2 

[mim2(f - e^i-^^)2e2^-^ + m2m3(l - e'^^^^'^fe^''^ + TOim3(l - e^i-^^)2e2^^] 

- [mim2m3(l - e^i-^^)2(l - e^^-^^)2e2^^] z^, (2.7c) 

The t evolution of A, B, C can easily be inferred from the second equation 
of the Lax pair (|2.ip . One obtains: 

B " 
i = 0, B = 2AM+, C = -BM+, where M+ = ^ TOic'^' . (2.8) 

We therefore see that the asymptotic conditions p.3p can be implemented by 
requiring A{z) = 0, and that this condition is preserved under the time flow 
(|2.8p . implying that the peakon equations with 6 = 3 describe an isospectral 
deformation of the boundary value problem (12. Sp . This boundary value problem 
can be best studied with the help of two rational functions 

Definition 2.1. (Weyl functions) u;(z) = Ci^) = 

In this paper we will only use U!{z). For the case of pure peakons, nii > 0, 
it was proved in [53] that uj{z) is a Stieltjes transform of a measure, which 
subsequently played a major role in the solution of the inverse problem. If, 
however, m is a signed measure then uj{z) has a more complicated structure 
because the spectrum is, in general, not simple or even real. Yet, nontrivial 
information about the dynamics of peakons can be extracted from uj{z) without 
knowing its precise pole structure. To this end let us establish a simple lemma 
which follows trivially from equations (j2.8l) and the definition of uj{z). 
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Lemma 2.2. 



. , , Lu(z) M+ 

u;(z) = -^ + — (2.9) 



From explicit formulas A{{)) = 1, B{0) = 0, thus implying that is a remov- 
able singular point of ^(2). Moreover, knowing the evolution of ^(2) we can 
readily establish the time evolution of the data involved in its partial fraction 
decomposition. 
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Theorem 2.3. Suppose the partial fraction decomposition of (jj{z) is given: 

% hf\t) 



J k=i ^ 

where dj is the algebraic degeneracy of the eigenvalue Xj . Then 

bf\t)^pf\t)e^, (2.10) 
where p'"^ {t) is a polynomial in t of degree dj ~ k or lower, and 

Y^hf\t)^M+. (2.11) 

3 

Proof. Combining the partial fraction decomposition with (j2.9l) one gets 

j k=l ^ ^' 

By Residue Theorem, we have = Res oo) +Res o) +Ej Res (^^, Xj 

where 

Res (^^,oo) ^ 0, Res (i^,o) = Res (^,A,) = ^^(t), 

which proves (|2.1ip . 

By the formulas for the coefficients in the Laurant series of equation (|2.12l) 
we obtain 



^ ^ ' ^ {s--k)\Az^-^ \ z 



s—k 



, s=k j 



In particular, we have b'"J^^\t) — -^^-y-^, hence b^^^\t) = h^^^\Q)e^j . Proceed- 
ing by induction we obtain 

s—k 3 s—k+1 j ^ 

therefore bf\t) = (bf\o) + Jo P^''^^HT)dT^ =^ e^pf\t). Finally, since 

^(fc+i) j-^) jg polynomial of degree rf^ — — 1 or lower, p^'^-' (i) is a polynomial 
of degree dj — fc or lower, which leads to p.lOp . □ 
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Lemma 2.4. Let Xn be the position of the n-th mass. Then 

e^-=Y,bf\ (2.13) 

3 

Proof. By Residue Theorem = Res (a;(z), oo) + Res (w(2), Aj). Thus 

^.jb^'p = — Res cx)). With the help of exphcit formulas (12.51) and the 

definition of lo{z) we obtain 

Res (aj(z), oo) — — lim zuj{z) — lini ^ \ — — e^", 
2— >oo z— >oo 2Ayz) 

which proves the conclusion. □ 

Corollary 2.5. The nth mass cannot escape to +oo in finite real time. 

Proof. Indeed, from Theorem 12.31 and the lemma above we see that e^" has at 
most an exponential growth, hence it is bounded for finite real time. □ 

To deal with the behaviour of Xi we will use a slightly modified spectral 
problem which, in principle, amounts to "sweeping" the masses in the opposite 
direction. To this end we consider the adjoint Lax pair: 

(a, - 9,x,)* = -zm^-, ^t = [-z-\l-dl)+u,-ud^]^. (2.14) 

Remark 2.6. The only difference between equations p.ip and (|2.14|) is the 
sign of z which has no effect on the compatibility conditions; hence the adjoint 
Lax pair gives the same compatibility condition — the DP equation. 

We choose a different set of asymptotic conditions, namely 

^ ^ e~^, as X — > +00, ^ is bounded as x ^ — oo. (2-15) 

For X < xi 

^{x)=A{z)e-'' + B{z) + C{z)e''. (2.16) 

Hence the adjoint spectral problem is given by A{z) = 0. Likewise, one readily 
checks that the time flow given by the second equation in (|2.14|) yields: 

i = 0, B = -— + 2AM_, C = BM_, where M_ = V m^e'''\ (2.17) 
z ^—^ 

We conclude that the adjoint boundary value problem (|2.15l) is also isospec- 
tral under the DP flow. In fact, the spectral problems (12. 3p and (|2.15p have 
identical spectra. To demonstrate that we establish first an elementary lemma. 

Lemma 2.7. // '^(x) is the solution to the x- equation in the boundary value 
problem (j2.1l) with m{x) — X^ILi "^^'^^^i' then '^(—x) is the solution to the x- 
equation in (|2.14[) with ifi{x) — "Y^^^i nT'iS-xi and boundary conditions (I2.15p . 
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Proof. Since ^(x) is the solution to (|2.ip . the boundary conditions 

^'(— x) e~^, as a; — > +00, 'i'{—x) is bounded as x — > — cxo 
hold. Moreover, we have 

{da: - dxxx)^{-x) = -{^ x{-x) - ^ xxx{-x)) = - Zm{- x)^ {- x) . 

Notice that a^) = ^-xi{^\ hence m{—x) — m{x) and the conclusion holds. 

□ 

Denote m = (toi, . • . , mn),x_ = (xi, . . . , a;„) for short, and set rn^,x^ to be 
the vector with the reversed order of its entries, that is rn^ = {nrim . . . , mi) etc. 
Employing the same convention as in equation (|2.4p . but this time for '5, we 
obtain the following analogue of equation (12. 5p . 

Theorem 2.8. Let Ak{z) = Ak{z;rn,x), 1 < k < n. Then 



Ck(z) J \0 J P=i 




0^ I I -2e-^'i 

-2xi-^ 



e 



where (^^p"') is the set of all p- element subsets I = {ii < • • • < ip} of {1, . . . ,n}. 
In particular, when k = n, A{z) — A{z). 

Proof. By lemma [2771 ^(x) in the asymptotic region x — > —00 can be expressed 
as 

^(x;m,a;) — A{z;ni,x)e'^^ + B{z;rn,x) + C{z;m,x)e^ 
=5'(-a;;rri^, -xj) = A{z;mJ', -x^)e"^ + B{z;rrf , -x^) + C{z;rn'', -x^)e^' 

which leads to A{z;m,x) — A{z;m'^ , —x'^), B{z;m,x) — B{z]frif , —xj'), and 
C{z;rn,x) = C{z]rnf , —xf). The conclusion then directly follows from the for- 
mulas (|2.5p along with an elementary observation that the permutation t is a 
bijection on the ordered p-tuples, which for any fixed p-tuple maps the last ele- 
ment e^*p in the original sum into the first element 6"^"+"^"'^ of the new p-tuple. 
After a simple change of index the main claim is proven. As to ^„(z), which 
corresponds to the first line in the formula, we observe that this polynomial is 
invariant under the transformation m 1-^ mf ,x_ ^ —xf. □ 

We can thus define the adjoint Weyl function oj(z) = — tt^tt^, and use 
equations (I2.17P to determine the time fiow of ui. An easy computation gives: 

Lemma 2.9. _ 

= ^ . (2.18) 

z z 



Consequently, we obtain an analogue of Theorem 12.3 
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Theorem 2.10. Suppose the partial fraction decomposition of uj{z) is given: 



J k=i ^ 3' 

where dj is the algebraic degeneracy of the eigenvalue Xj . Then 





x(fc) 



^'^'^ -Pf\t)e~^ (2.19) 



where (t) is a polynomial in t of degree dj — k or lower, and 

^^f^=-M_. (2.20) 



With the help of Theorem 12. 81 it is now straightforward to establish a coun- 
terpart of Lemma 12.41 

Lemma 2.11. Let xi he the position of the first mass. Then 

e-^^ ^J2f\t). (2.21) 

i 

Likewise, an analogue of Corollary [53] is immediate. 

Corollary 2.12. The first mass cannot escape to —oo in finite real time. 

Example 2.13. Case n = 3. In this case we are only dealing with simple and 
quadratic roots, since the triple roots cannot occur as will be proved in Section 
4. The formulas for e^^^*^ and e^^i^^*^ read: 

b[^'' (O)e^ + (O)e^ + 5^^^ (O)e^ , simple roots, 
b[^'> (O)e^ + (4'' (0) - ^^^|^)e^ , quadratic roots. 
4^^(0)6"^^ +4^^0)6^^^ +4^^(0)e~'^, simple roots, 
b[^\0)e^^ + (4^^(0) + ^ quadratic roots. 

The spectral problem and its adjoint are clearly related and we turn now 
to establishing a relation between them. To this end we study the coefficients 
occurring in the eigenfunctions of the spectral problem (|2.6p and (|2.16p . 



Lemma 2.14. 



2A{z)C{~z) + 2A{~z)C{z) - B[z)B{-z) = 0, 
2A{z)C{-z) + 2A{-z)C{z) - B{z)B{-z) = 0. 
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Proof. It suffices to write {D — D^)ip{x;z) = zmip{x; z), {D — D^)ip{x; X) = 
Xmipix; A) and obtain from it the identity: 

D{^{x; z)ip{x; A) — B{tp{x; z),ip{x; A))) = (A + z)mip{x; z)ip{x; A), 

wliere B{f, g) = f" g — f'g' + fg". Finally, if one sets z + X — and evaluates 
the above expression at cc — — oo and x ^ oo one obtains the first claim. The 
proof of the second identity is analogous. □ 

We will briefly study the symmetry responsible for the connection between 
the boundary value problem (|2.3p and its adjoint (|2.15p . To this end we recall 
the transition matrix S{z) introduced in [23] 









'Ak-i' 


S{z) = Sniz)Sn-i{z) ■ ■■ Si{z), whcrc 


Bk 


= Sk{z) 


Bk-i 




Ck 







where Aq = 1, Bq = Cq = 0. An explicit form of Sk{z) is easy to compute: 



Sk{z) = I - znik 

Define now 
Definition 2.15. J 



(2.22) 



1 
0-2 

1 



We can define the loop group of continuous maps G : R ^ SL(3,R); 
clearly S{z) G G. Moreover, if we introduce involution: r : G — > G,g{z) — !> 
J((7^^(— z))^ J^^, then Sk{z) E Gr = {g = T{g)}, a subgroup fixed by r. Hence 

Lemma 2.16. S{z) e Gr- 



Let us denote the canonical basis ei = 



'I 




"0" 




"0" 





,62 = 


1 


,63 















1 



by the 



shorthand notation |2), |3). Then equations (|2^ . (|2?T6l) can be written 



A{z) 
B{z) 
C{z) 



= S{z)\l) 



■C{z) 
B{z) 
Mz). 



= 5(z)|3). 



where 



S{z) = S^\~z) ■ • • S-\-z) - S{-z)-\ 
Moreover, since Sk{z) G G,-, 

S{z) = JS'^{z)J-^. 



(2.23) 

(2.24) 
(2.25) 



This is a fundamental relation which allows one to relate the spectral data for 
the boundary problem and its adjoint. 
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Theorem 2.17. 

J. = = -25(z)i2, C'(z) = 5(z)i3. 

2. Suppose Xi is a root of A{z) — then 

B{-\i) = C{\i)B{\,). (2.26) 

3. Suppose B{—\i) ^ then 

2A{~\) = B{\i)B{X{). (2.27) 

Proof. By definition A{z) = {i\J {z)J-^\i) = (l|5'^(z)|l} = S'(2)ii. Like- 
wise, B{z) = (2|J5'^(z)J-i|3) = -2(2|5^(z)|l) = -25(z)i2 and C{z) = 
l^\J S'^ {z)J^^\Z) = 5(z)i3. The second item is proved by making use of the 
involution r. On one hand B{z) = (2|S'(z)|l), on the other, since S{z) G Gt, 

B{z) = -2{3\S-\z)\2) - 2 ^J^jj" . Finally, since 5(z)n = A{z), 

evaluating the determinant at the (minus) root Xi of yl(z) we obtain B{'-Xi) = 
— 25'(Ai)i25'(Ai)3i = B{Xi)C{Xi), in view of the statement from item (1). Fi- 
nally, to prove item (3), we set z = Ai in the statement of Lemma [2.141 to get 
2A(-A,)C(A,) = B{Xi)B{-Xi). Upon multiplying equation (12:26)) by 2A(-Ai) 
and eliminating the term involving 2A{—Xi)C{Xi) we obtain 2A{—Xi)B{—Xi) = 
B{-Xi)B{Xi}Bi{Xi), resulting in equation (P??fl) . □ 

Consider now the Weyl function u){z) and its adjoint u}{z) in the case of 

simple spectrum, i.e. 

B{z) ^ bi 1 -/ N B{z) ^ &j 

uj(z) ^ -- — — — = > — and lu(z) = = > — . 

^ ' 2zA{z) ^ z - A, ^ ' 2zA{z) z - X, 

Theorem 2.18. Suppose the spectral problem A{z) = has only simple roots 
Xi and there are no anti-resonances {Xi + Xj ^ 0). Then 



^^^^ = n71^^- (2-28) 



Proof. Under the assumption of simple spectrum: 

~ _ B{Xi)B{X,) 
- 4AKA'(A0)2 

which simplifies, after using equation p.27p . to 



bibi = 



2xnA'{x,w 2n,^.(i-^ 



which implies the claim if one observes that the term with j = i appearing in the 
numerator contributes the factor of 2 canceling the one from the numerator. □ 
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Remark 2.19. This beautiful identity generalizes the one known from the 
ordinary string problem |25| which in our notation reads: 

Remark 2.20. The presence of anti-resonances (Ai + Xj = 0) is characteristic 
of the DP equation as can be seen, for example, from explicit solutions. 

3 Three multipeakons 

In this section we apply the methods developed in Section [2] to study three 
multipeakons, with emphasis on the behaviour of solutions at the time of blow- 
up. As before we use the multipeakon ansatz (|1.4p 

3 

= ^m,(t)e-l^-^'(*)l (3.1) 
1=1 

where a;i(0) < 2:2(0) < 2:3(0), and we no longer assume that mi{t) are all 
positive. In spite of that we will refer to m^s as masses to emphasize their 
roles in the spectral problem. We will need a bit of terminology regarding the 
phenomenon of breaking. Since we will be analyzing a system of ODEs obtained 
from a restriction of equation (|1.7|) we will say that at some time to a collision 
occured if for some i ^ j, Xi{to) = Xj{tQ). In the case of the CH equation the 
presence of a collision is tantamount to a wave breaking ( [14 ) but the solution 
can be continued with the preservation of the Sobolev i/^(R) norm beyond the 
collision time. This is not the case for the DP equation as was anticipated by 
Lundmark in [TT] for the case of the peakon-antipeakon pair. We confirm his 
assertion that the shockpeakons are created by proving that m — u — u^x tends 
to the shockpeakon data in the distribution topology at the collision time (see 
Theorem Eini)- 

We start by setting 6 = 3 and n = 3 in the multipeakon equation (|1.7p . 
which leads to the following ODEs in the sector X = {x G R"^ |a;i < < x^}: 



±1 


= mi -f m2e'^^^'^^ -|- m3e'^^~'^^. 


(3.2a) 


X2 ■ 


= mie"^^""^^ -1- m2 -1- m3e'^^~'^^. 


(3.2b) 


xz 


- mie^^^^^ + 17126^^^^^ + m3, 


(3.2c) 


mi 


= 2mi(-m2e^i~^^ - mae^i"^^'). 


(3.2d) 


m2 


= 2m2(mie'^i^'^^ - m^e^^''^^), 


(3.2e) 


7713 


= 2m3(mie^i""^3 + m2e="=^-^='). 


(3.2f) 



This system of equations has the following obvious symmetry. 

Lemma 3.1. Suppose {xi{t),X2(t),X3(t), 1711(1), m2(t),'m^{t)} is a solution of 

equations (13. 2p at time t with the initial condition {xi(0), 2:2(0), 2:3(0), mi(0), m2(0), m3(0)}. 
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Then {xi{t),X2{t),X3{t), —mi{t), —m2(t), —msit)} is the solution at time —t 
with the initial condition {xi(0), 2:2(0), 2:3(0), —mi (0), —m2(0), —771,3(0)}. 

Remark 3.2. In short, the lemma above means that t n> —t,mi 1-^ —rrii is a 
symmetry of equations p.2p which preserves the sector X. 

Another very useful property of equations (13.21) is the existence of three 
constants of motion. Indeed, we recall that the polynomial A{z) introduced in 
(|2.7I) is time invariant. Writing 

A{z) = 1 - Miz + Khz"^ - M^z^ (3.3) 

we obtain the following lemma. 

Lemma 3.3. Mi, M2, M3, given by: 

Ml =7771+7772 + 7773, 

M2 = mi7772(l - e"^!""^^)^ + 77727773(1 - e^"^-^^')^ + 777377n(l - e"^l""^^')^ 
M3 = 777177727773(1 - e'''-'=^f{l ~ e='''-='^f, 

are constants of motion of the system of equations p. 21) . 

These constants will be one of our basic tools for studying collisions. We ob- 
serve that, geometrically speaking, a collision occurs if the solution approaches 
the boundary of X in finite time. This is the only singular behaviour of the 
system (|3.2I) happening in the coordinate space since Corollaries 12.51 and 12.121 
exclude an escape scenario in finite real time. However, the shape of the con- 
stants of motion shows that at a collision at least two masses diverge, which will 
be proved in Corollarv l3.16l 

We begin now our study of the dynamics of three multipeakons in a vicinity 
of the collision by first concentrating on the particles with labels 1 and 3. Lemma 
12. 4[ in particular equation (|2.13p , gives us explicit form of X3 : 

gX3(t) (3.4) 
3 

Thus we obtain: 

Lemma 3.4. LetTi, T2, he the largest negative, respectively the smallest positive 
root of'^jb^^\t) =0 (if Ti or T2 does not exist we set Ti = — cxD,r2 = +00 
respectively). Then e^^*-*-* is real analytic for Ti <t<T2. Moreover, if either Ti 
or T2 are finite then there must be a collision at some prior time Ti < tc < T2. 

Remark 3.5. For positive t "prior" has the usual meaning (positive orienta- 
tion). For negative t the orientation is from to —00. 

Proof. The last statement follows from Corollary 12.121 since finite Ti or T2 
means that 2:3 escaped to —00 which cannot happen in finite time unless there 
is a collision at an earlier time. □ 
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Remark 3.6. Since b^^Ht) is an exponential function of t, the right hand side 
of p.4p is well define for any real t. However the left hand side of p.4p only 
make sense when t lies in the existence interval of the ODE system p.2|) . 



Furthermore, combining equations for ±3 and (see equations (j3.2p ') yields: 

7713 = 2m3(i3 - TO3). (3.5) 

We remark that this is a Bernoulli type equation which can be easily solved 
once x^lt) is known. 

Lemma 3.7. Suppose X3{t) is known. Then 

(3.6) 



1 r„2a;3(0) ft 



^ e -2a;3(t) 

TO3(i) Lm3(0) JO 

An analogous argument works for xi. Indeed, by equation (|2.21|) . we have 



This prompts an analogous statement to Lemma 

Lemma 3.8. Let Ti, T2, he the largest negative, respectively the smallest pos- 
itive, root of ^jU^^t) = 0. Then e"^^'*-* is real analytic for Ti < t < T2. 
Moreover, if either Ti or T2 are finite then there must he a collision at some 
prior time Ti < tc < T2. 

We see that we can now narrow down the time of a collision. Let us denote 
hy A = {Ti,T2) n (Ti,T2). We summarize analytic properties of e^^ and e~^^. 

Lemma 3.9. The functions e'^^^*'^ ande"^^^*"^ are real analytic on A. A collision 
can only occur at a time tc if tc ^ A . In particular, both functions are analytic 
at the time of collision. 

Once again, if we know xi then we can determine mi. 
Indeed, equations p.2p imply another Bernoulli equation: 



mi — —2mi{xi — mi), (3.7) 



whose solution reads 
Lemma 3.10. 

1 



-2xi{0) ft 



2 / e-2-iMdr 



mi{t) I mi(0) JO 

We can now summarize analytic properties of mi,m2 and m3. 
Lemma 3.11. 



(3.8) 
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(1) j^^^ji^ and j^^f^f.-^ are real analytic on A. 

(2) m2{t) and e^^^*-* are real meromorphic functions on A. 

(3) A collision occurs iff there exists tc ^ A such that either 
—TT-x = 0- 

ni3{tc) 

(4) Suppose ^^^^^ J — 0. Then in a neighborhood of tc 

-J— 2it~t,)+oiit^t,r). 

mi{t) 

0. Then in a neighborhood of tc 

= 2{t-tc) + 0{{t-tc)^). 



or 



(5) Suppos 



1 



m^it) 

Proof. The analytic properties of — and — are directly derived from the 
analytic properties of e^^*^*^ and e^^^*^ and Lemmas 13. 10) [3771 respectively. 

To prove (2) we note that using Mi we can write m2 = Mi — (mi + ma). 
Likewise e^^ can be easily computed from M+ = nije^' = ^e^^ by algebraic 
operations on analytic functions. 

To see (3), since is a constant of motion, we observe that a collision 
occurs when at least two of the masses diverge, hence either mi or ms have to 
diverge at a collision. 

and 



To prove (4) and (5) we only need to calculate the derivatives of 



21 (t) 



^^(^^^ at tc- Since xi,X3 are analytic at tc, direct computation from (j3.5p and 
p.7p shows that 



1 



dt ms (t) 
d 1 



- 2- 



2X3 

ma 



dt mi{t) 
Therefore (3) and (4) hold. 



2ii 
mi 



- 2 



= 2, 
= -2. 



□ 



Now we can obtain the behaviour of masses before a collision. As a general 
comment, we observe that for any initial data in X and arbitrary mi,m2,rn3 
the solution is unique. From this point onwards we assume mi(0) ^ 0. Thus 
M3 is nonzero. 



Lemma 3.12. None of the masses mi can become zero before a collision. 



15 



Proof. In order to derive a contradiction, we suppose that one of the masses 
becomes zero at to. Since all three constants of motion given by Lemma [3731 are 
symmetric with respect to permutations of masses, we can assume, without a 
loss of generality, that mi{to) = 0. Since M3 ^ 0, raivn-i diverges at ^o- Then 
jj^^^ converges to zero, while the corresponding right hand side converges to 
(1 — e^^~^^)^ 7^ 0, thus a contradiction. □ 

Corollary 3.13. The masses rrii cannot change their signs before a collision. 

Corollary 3.14. None of the masses mi will diverge to ±00 before a collision. 

Proof. Since none of can become zero by Lemma 13.121 is nonzero, < 

M3 ^ J M: 



< 00 and — — = (1 — e^i ^'^Y{1 — e^^ ^'^Y is continuous before a 
collision, hence the claim follows. □ 

Combined with the analytic property, there are two corollaries worth men- 
tioning. 

Corollary 3.15 (Absence of triple collisions). There are no triple collisions, 
that is, there is no time at which xi — X2 = X3. 

Proof. Suppose xi(tc) = 2^2 (^c) = x^itc). Then the leading contribution to M3 
coming from the term (1 — e^^~^^Y{l — 6^2-2:3 -^2 (^((i _ tcY) which forces 
TOim2m3 to behave like 0{jj—^^jf). If only mi,m2 diverge then mi diverges as 
2(^1^ _f) by Lemma 13.111 and m2 would have to diverge as 0{ (j-J^^yi ) violating 
conservation of Mi . Similar argument excludes divergence of m2 , TO3 . The last 
case is that mi, m2, TO3 diverge, but then in view of Lemma l3.11l m9 would have 
to diverge as again violating conservation of Mi. □ 

Corollary 3.16. At the point of a collision masses diverge in pairs and the 
only admissible pairs are {mi, 7712} and {m2,TO3}. 

Proof. In view of the behaviour of mi and 7713 at the collision, m2 must be 
regular to preserve Afi if mi and TO3 diverge. Thus mi, 777,2,7713 cannot all 
diverge. To eliminate the mi,m3 pair we consider — limt_5.t^ nnm-j ^ ^ 
^xi(t^)-x3(t^)y2 ^ Q |-,y ^}jg absence of triple collisions, hence a contradiction. □ 

Theorem 3.17 (Shockpeakon creation). // mj collides with m^+i at tc > 0, 
then 

lim {mj{t)S{x — Xj{t)) + mj+i{t)5{x — Xj+i{t))) 
lim (mj + 77ij+i) ) S{x - x{tc)) + - ( lim {u{xj{t),t) - u{xj+i{t),t)) \ 5'{x - x{tc)), 
where the limit is in the sense of S>'{M}j. 
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Proof. For arbitrary Lp{x) 6 ^(M), 

{mj{t)S{x-Xj{t))+mj+i{t)5{x-Xj+i{t)),ip{x)) = mj{t)ip{xj{t))+mj+i{t)ip{xj+i{t)). 
Whenever j — 1 or 2, we can always write 

1 „ . 1 



+ Co + 0{t-tc), mj+i = — — — + Co + 0(t - tc) 



around t^. Hence, 

lim {■mj{t)S{x — Xj{t)) + ■mj^i{t)S{x — Xj^i{t)), (p{x)) 

lim_(TOj +mj+i) ) tf{x{tc)) - \[ lim_(xj -i^+i) ) ip'{x{tc)) 

lim {nij +mj+i) ) (p{x{tc)) - ^ ( lim {u{xj{t),t) - u{x j+i{t) ,t)) \ ip'{x{tc)), 

where in the last step we have used equation (|1.7p for & = 3. The claim follows 
now easily from the definitions of distributions S and S' . □ 

Remark 3.18. Shockpeakon creation described by Theorem 13. 1 71 confirms the 
scenario that at the collision the colliding peakon-antipeakon pair creates the 
shock (the 6' contribution above) and the peakon or antipeakon contribution 
(the S contribution) thus giving the overall collision data of two peakons/antipeakons 
and a shock. This has been previously verified for the case n = 2 in [TT]. 



4 Three multipeakons; spectral properties 

This section addresses basic questions related to the spectral characterization 
of the peakon dynamics (|3.2I) . 

Lemma 4.1. Let denote the number of positive masses and n+ be the 
number of eigenvalues of the spectral problem A{z) — \^ Miz+ M2z'^ ~ M^z'^ — 
which have strictly positive real parts. Then 

N+ = n+. 

Proof. The statement holds true if = 3 by results in [23]; in that case the 
spectrum is positive and simple. Since rrii H> —mi,Xi i— ^ — Ai is a symmetry 
of the eigenvalue problem it suffices to analyze only the case with two positive 
masses, that is — 2. 

Then ^ 

A1A2A3 
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To prove the claim we have to exclude that three eigenvalues have strictly neg- 
ative real parts (recalling that complex roots must occur in conjugate pairs) or 
that there is one negative and two purely imaginary conjugate roots. In either 
case Mi = Y- + T- + T-<0, and M2 = t^- + tV + t\- > 0. So we assume 

Ai A2 A3 ^ A1A2 A2A3 A3A1 

that Ml < and M2 > in order to derive a contradiction. 

Case 1. If mi < 0, then 1712,1713 > and < m2 + < —mi. Therefore 
-mim3(l - e^i-^^)2 > mamall - e^=-^^)^ 

which implies 

M2 = mim2(l - e^i-^=)2 + m2m3(l - 6^="=^=^)^ + miTO3(l - e^'^-'^' f < 

and thus leads to a contradiction. 
Case 2. If ma < 0, it is similar to Case 1. 
Case 3. If m2 < 0, then mi, ma > and < mi + ms < —m2- 
Denote m = mi + ms, then 

mim2(l-e^i-^=)2+m2m3(l-e^^-^=')2 < -m[mx{l-e''-''^f+m3{l-e'^-'=^f]. 

Set a = e^l-^^/3 = e^^-^^ml = 6im,m3 = (1 - 6»)m, and 

f(6) = 6{\ - af + (1 - e){\ - - e{i - e){\ - aPf. 

f{9) is a quadratic function with respect to with the discriminant 

A = ((1 - af - (1 - (3f - (1 - - 4(1 - a/Sfil - Hf 

= -(1 - a)2(l - /3)2(l + a){l + /3)(3 - a - /3 - a/3) < 0, 

which leads to f{9) > 0. Therefore 

M2 < -m[mi(l - e^i-^=)2 + m3(l - e^^"^^)^] + m3mi(l - e^^"^^)^ 
= -m'^{e{l - af + (1 - e){l - - 6(1 - e){l - apf) = -rn^fie) < 0, 

hence a contradiction. □ 

Clearly, by reflection symmetry, we obtain 

Corollary 4.2. Let N~ denote the number of negative masses and n~ be the 
number of eigenvalues of the spectral problem A{z) = l—MiZ+M2Z^—M3Z^ = 
which have strictly negative real parts. Then 

N- =n-. 

Another useful corollary is that there are no eigenvalues on the line Re 2; = 0. 

Corollary 4.3. None of the eigenvalues of the spectral problem A{z) = 1 — 
Miz + M2Z^ — Mzz^ = is purely imaginary. 
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Corollary 4.4. The spectral problem for n — 3 can never have triple roots. 

Proof. Suppose, without loss of generality, that the spectral problem has triple 
positive roots. Then all the masses are positive, i.e. the peakons case. However 
the eigenvalues for the peakons are simple, hence a contradiction. □ 

Remark 4.5. Figure[T]on pagefTOlillustrates how the eigenvalues are distributed 
for the mass signature mi(0) > 0, 77x2(0) < 0,7713(0) > 0, abbreviated 

(H h). The graph depicts 75 x 75 triples of eigenvalues for different values 

of masses within that configuration. The actual input data is mi — 1.2 + 
0.02j, 7772 = -5 - O.Olfc, 7773 = 4, Xi = -0.2, X2 = 0, X3 = 0.1, I < j,k < 75. 
Observe that indeed the line Re A = contains no eigenvalues. 



6 




_ei 1 1 1 1 1 

-15 -10 -5 5 10 



Figure 1: A portrait of eigenvalue distribution for the mass signature (H h) 
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5 Three multipeakons; classification 



The goal of this section is to understand the impact of the configuration of signs 
of masses on the occurrence of coUisions. 

One can classify the system of three multipeakons by the signs of the initial 
values of miS: 

(i)TOi(O) > 0,TO2(0) > 0,m3(0) > 0; (ii)mi(O) > 0,TO2(0) > 0,TO3(0) < 0; 
(iii)mi(O) > 0,m2(0) < 0, 7713(0) > 0; (iv) 777i(0) > 0,7773(0) < 0,7773(0) < 0; 

(v) 777l(0) < 0,7772(0) > 0,7773(0) > 0; (vi) 777l(0) < 0,7772(0) > 0,7773(0) < 0; 
(Vii) 777l(0) < 0,7^2(0) < 0,7713(0) > 0; (Viii) 777l(0) < 0,7r72(0) < 0,7773(0) < 0. 

The first and last cases are pure peakon and antipeakon, which are already well- 
known. In view of Lemma [3.11 the symmetry 777^ — >■ —rrii^t — >■ —t reduces the 
eight cases to four cases. To gain some clarity we will supplement a reference to 
any of the cases from the list above by an ordered collection of signs, i.e. case 
(i) is equivalent to (+ + +), case (ii) to (+ H — ) etc. and we will refer to a given 
mass signature as a 7i7ass signature. 

Theorem 5.1. If the mass signature is: 

(1) ( h) or { — h +), then no collisions will occur for positive times, 

(2) (+ H — ) or (H ), then no collision will occur for negative times. 

Furthermore, if the eigenvalues are not in anti-resonance, a collision will always 
happen at some finite time tc 

(1) if the mass signature is ( h) or ( — h+), then the collision will happen 

at a negative time, 

(2) if the mass signature is (+ H — ) or (H ), then the collision will happen 

at a positive time. 

(3) if the mass signature is ( — \ — ) or (H h), then the collision will happen 

at both a finite positive and a finite negative time. 

Proof. First, we show that in the case of item (1) no collisions occur in positive 
time. Indeed, by examining the formulas (|3.6p and (|3.8|) we see that the respec- 
tive right hand sides can not be equal for t > 0. The same argument works 
for item (2) and negative times. 

Let Ai , A2, A3 be the eigenvalues of the system which we can order as Re < 
Re^ < Re^. Since the masses have different signs, by Lemma l4.ll we have 
Re^ < < Re^. 

Case 1 : The eigenvalues are simple. Since the eigenvalues are not in anti- 
resonance, all the residues 6i's and hiS are nonzero according to (|2.28p . Hence, 
according to (I2.13P and (I2.2ip . there exists at least one increasing and one 
decreasing exponential function in the expansions of e^^ (respectively in the 
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expansions of e~^^). Moreover, since we are squaring e^^, e~'^^ respectively, the 
coefficient of the leading exponential will be strictly positive if the spectrum 
is real, or strictly positive except for a set of measure zero if the spectrum is 
degenerate or complex. That is to say, both integrals 

t nt 

e^^^^^Mrand / e^^^^^^^Mr 

JO 

will diverge to icxi as < — > ±00. Hence, there exists a positive (respectively 
negative) time such that ^ = or ^ = whenever 7713(0) < or mi(0) > 
(respectively 7713(0) > or 7Mi(0) < 0). This proves the claim in view of lemma 

Case 2: There is a double root, i.e. Ai = A2 ^ A3, then Ai will not be the 
double root of B{z) (respectively B{z)), therefore at least one of b^i\o) and 
^^^■'(O) (respectively b^j^\o) and ^^^■'(O)) is nonzero. This also implies that there 
exists at least one increasing and one decreasing exponential function in the 
expansions of both e^^ and e~^^ . To show that Ai will not be the double root of 
B{z), we only need to note that Ai must be the double root of A{—z)C{z) if Ai 
is the double root for both A{z) and B{z) according to Lemma [2. 141 However, 
C(Ai) must be nonzero, otherwise vE'(z) is identically equal to which leads to a 
contradiction. Therefore — Ai must be the double root of A{z), which also leads 
to a contradiction. □ 

It is immediate from the above theorem that 

Corollary 5.2. 

(1) For cases (v){ h +) and (vii){ h) there exists a unique solution to 

the ODEs p.2p for all positive t. 

(2) For cases (ii){+^ — ) and (iv){^ ) there exists a unique solution to the 

ODEs (j3.2p for all negative t. 

The global existence (in one time direction) has an interesting impact on the 
spectrum of the boundary value problem. 

Theorem 5.3. The eigenvalues of the spectral problem of cases (ii,iv,v,vii) are 
real, simple, nonzero, and are equal to the inverses of asymptotic values of 
masses. 

Proof. We give a complete proof for the case (ii)(+ H — ). First, we note that 

m3(i) = m3(0)exp 1^-2 y"[77ii(s)e^i('')-"=^('') + m2(s)e^^('')-"=^('')]ds^ ^Q. 

For negative t, recalling that mi, 777.2 will remain positive, we obtain 

|m3(t)| < |m3(0)|. 
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Since mi < mi +m2 < Mi + \m3(t)\ we derive an upper bound on mi, and thus 
on 7712, namely 

m,{t) < Ml + \m3{0)\ , 1^1,2. 

We observe that < Mi + |m3(0)|, otherwise mi — m2 — for all times. 
Claim I 

^1 ' iM ^ <m,(t) <Mi + |m3(0)|, t<0. 

[Ml + |m3(0)|) |m3(0)| 

Proof. (Claim I) We only need to prove the lower bound. To this end we esti- 
mate: 

IM3I < mi{t)m2{t) \m3{t)\ < mi{t)m2{t) |m3(0)| , 
and use the upper bound above on one of the factors mi or m2. □ 

Replacing in the estimate for IM3I both factors mi and m2 with their upper 
bounds we extend the claim to the bound on m3. 
Claim II 

Claim III 

Xi{t) — Xj{t) —00, (1 < i < J < 3), 

when t — ^ —00. 

Proof. (Claim III) The following estimate holds: 
IM3I 



(Afi + |m3(0)|)2|m3(0)| 



< 







exp(-2 / [mi(s)e^i(')-^-'(") +m2(s)e^^('')-^-'(")]ds) 

and thus 



[mi(s)e^i(")-^^(") +m2(s)e^^(")-^-^(")]ds < 00, 



which, in view of the boundedness of mi and m2, implies 




e^^'f^^-^^^'^Ms < +CX), i = l,2. 



In addition, direct estimates on equations p.2p using the bounds on mi, m2, m3, 
show that the velocities are bounded, which means the derivative of the inte- 
grand e^i(«)-^3(s) ^-j _ ]^^2) is bounded. Therefore 

lim e^'(")-^^(*) =0, (i = l,2). 
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which is equivalent to 

hm Xi{t) — X3{t) = — oo. (i ~ 1, 2). 

Now we turn to the ODE for mi 

mi{t) = mi(0) exp I -2 / [m2is)e''''^''^~'='^'^ + m3{s)e''''^''^-'='^'^]ds I > 0. 



Since nii is, for negative t, bounded from above and from below away from 0, 
the integral 





and J e^^*-*-* ^^(^^ds < +oo. Repeating verbatim the arguments from the pre- 

— oo 

vious case, we get 

lim xi{t) — X2{t) = — oo. 

□ 

Now, since the improper integrals appearing in the formulas for mi and ma 
exist, we can take the limit t — >■ — oo. Let us denote those limits by mi(— oo), 
m3(— oo) respectively. Using Mi we conclude that m2 also has a limit, say, 
m2(— oo). The characteristic polynomial A{z) (see p.3p ) reads: 

A(z) = (l-f )(l-f )(l-f ) 

Ai A2 A3 

= (1 — mi(— oo)z)(l — m2(— oo)z)(l — ?7i3(— oo)z). 
Claim IV limt_^_oo a^i = linit^-oo X2 ^ - limf^_oo 2^3 = -00. 

Proof. (Claim IV) We prove first the claim for xi. The right hand side in the 
equation for ±1 (see p.2p ) reads mi + m2e^^^^^ + mze^^^^^ . When t —00 
the second and third terms go to 0. The limit mi (—00) > because of the lower 
bound on mi. Hence there exists a constant a > and another constant t* <Q 
such that 

0<Q;<ii, for aU t < i*. 

Integrating this inequality from (negative) t to t* we obtain xi{t) < at + C, 
where C is a constant, which proves that xi — t- —00 as t — > —00. The same 
argument works for X2- Since m3(— 00) < we get an opposite estimate for 0:3, 
namely, there exists a constant (3 < and another constant t** < such that 

is < /3 < 0, for all t < t** 

which, upon integration, yields l3t + D < 0:3 (t), forcing X3{t) — s> +00 when 
t -00. □ 
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Since 77i3(— oo) < 0, to prove simplicity, we need to show that the remaining 
two positive hmits mi (—00) and m2(— 00) are distinct. We adapt the proof of a 
similar statement in [53]. First, we observe that in view of (13. 2p m2 > 0, hence 
TO2(i) is increasing. The right hand side of the equation for mi reads (after 
dividing by 2mi) —m2e^^~^'^ — ■mze^^~^^ = e'^^~^'^{—ni2 — mze^^~^^). Since 
the second term goes to as i — —00 , and mi is bounded away from 0, we 
obtain that there exists t* <Q such that: 



mi < 0, for aU t <t* , 



So mi is decreasing for t sufficiently large and negative. Suppose, to derive a 
contradiction, that mi (—00) = m2(— cxd) then 

mi(t) - m2(t) < 0, foralH<r. 

On the other hand, 

M 

(xi ~ X2)dT ^ Xi{G) ~~ X2(Q) ~ {xi{t) " X2{t)) = J {mi ~ m2)d T + B (t) , 

where B{t) has a finite limit as t — > —00. The left hand side diverges to +00 
as t —00, while the right hand side can only diverge to —00 based on the 
inequality above. This contradiction shows that mi(— cxd) > m2(— 00). 

The proof for the other three cases is analogous. □ 

Remark 5.4. It is helpful to have an intuitive understanding of the above 
theorems. The emerging picture is this: if one has a swarm of peakons colliding 
with a swarm of antipeakons then the system essentially behaves as if it were a 
peakon-antipeakon pair. Thus the system has the following characteristics: 

1. peakons and antipeakons are asympotically (in an appropriate time direc- 
tion) free with asymptotic velocities xj = j- which are distinct 

2. peakons, antipeakons separate, that is Xi — Xj —00 for i < j. 

There are only two cases left: ( — I — ) and (H h). In view of the reflection 

symmetry (see Lemma [3. it suffices to analyze only one of them. We choose 

(- + -)• 

Theorem 5.5. Suppose the mass signature is ( 1 ), then there exists a 

positive time tc such that xi{tc) < X2{tc) — 2^3(^0) and a negative time t* such 
that xi{t*) = X2{t*) < X3{t*). 

Proof. We only need to prove that mi , m2 will never collide at a positive time. 
If not, assume that there exists a positive tc such that mi — )■ — cx),m2 — >■ +00 
and ma remains bounded when t tc- Then by equation (|3.2|) 

^(xi - X2) = (1 - e"i-"^)(mi - m2 - mge"^-"^) < 0, 
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when t is sufRciently close to tc- This contradicts xi — X2 — > 0. Likewise, for 
negative times, we need to ehminate a cohision of X2 and x^. However, 

^^{X2 - X3) = (1 - e"^^"^)(m2 - mg) - mi(e"i-"^ - e^^-"^) > 

for sufficiently close to the collision time < 0, again contradicting that X2 — 
as t ^ t*+. □ 

Below we put the results of our investigation in a form of Table [T] Speaking of 
the asymptotic behaviour we denote by AF the system which is asymptotically 
free in both time directions, if it is only in one, say in the direction of positive 
time, then we abbreviate it as AF-I-, etc. 



Table 1: Correlations of masses, spectrum and asymptotics 



Mass signature Spectrum 


Asymptotic behaviour 


Collisions 


+ + + + + + 


AF 


none 


+ + - + + - 


AF- 


< <c 


H h Ai < < Re A2 < Re As 


confined 


t* <0 <tc 


H h 


AF- 


< tc 


- + + - + + 


AF+ 


t* < 


- + - Re Ai < Re A2 < < A3 


confined 


tl<0<tc 


h - - + 


AF+ 


t*<0 




AF 


none 



We would like to conclude this section by discussing briefly the question of 
spectral data and in which sense the formulas obtained in J23j can be used to 
produce multipeakon solutions. We recall, in the notation of that paper, 

Xk'^log- — , ruk' = (fc = l,...,n), (5.1) 

vfc-i yVkWk-i 

where fc' = n + 1 — fc, and Uk and Vk are certain rational functions of spectral 
data 

{(A, 6) e R2" : < Ai < ■•• < A„, aU b, > O}. (5.2) 

(see formulas (2.44) and (2.45) in [23 for definitions), while Wj = UjVj - 
Uj+iVj-i. Clearly, for multipeakons, 

1. the spectrum is no longer positive, simple, or even real, 

2. the residues bj can be negative and, in general, complex, 

3. the anti-resonance condition Ai + Aj =0 renders the formulas not directly 
applicable 

4. once cannot extend the formulas beyond a collision point because there is 
no longer guarantee that Xi < Xj for i < j. 
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(a) DP mass signature (+H — ); 7712 collides with at tc = 0. 
Crossing of X2 with X3 




(b) CH mass signature (+H — ); m2 collides with ma at tc = 0. 
Collision of m2 with ms, followed by another collision of 1712 
with mi. No crossing. 



Figure 2: Comparison of DP and CH collisions 
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Example 5.6. Figure 2(a) on page illustrates how the formulas would work 

for the mass signature (+ H ). A point from the spectral set is chosen so 

that the collision occurs at tc = 0. The continuation beyond the collision point 
would force new ordering xi < X3 < X2 which means that this is not the original 
multipeakon problem given by equation (ll.Sp for & = 3, n = 3, even though the 
solution still satisfies equation p.2p , albeit in the wrong region. This should be 
contrasted with the behaviour of peakons at collision points for the CH equation 
as illustrated by figure [2(b)] on page The second particle bounces between 
77x3 and nil and, consequently, no change of ordering is required. 

Example 5.7. In this example we consider the mass signature (H h). Figure 

[3] on page E71 illustrates how the formulas (|5.ip would work for the mass signa- 
ture (H h). In particular, as predicted by Theorem 15.51 fusing the reflection 

symmetry), there are two collision points, one for negative time, one for positive. 




Figure 3: DP mass signature (H h); m2 collides with 7713 at < while 7712 

collides with mi for tc > 0. Confined state. 
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